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Multimodal signals facilitate communication with conspecifics during court- 
ship, but they can also alert eavesdropper predators. Hence, signallers face 
two pressures: enticing partners to mate and avoiding detection by enemies. 
Undefended organisms with limited escape abilities are expected to mini- 
mize predator recognition over mate attraction by limiting or modifying 
their signalling. Alternatively, organisms with anti-predator mechanisms 
such as aposematism (i.e. unprofitability signalled by waming cues) might 
elaborate mating signals as a consequence of reduced predation. We hypoth- 
esize that calls diversified in association with aposematism. To test this, we 
assembled a large acoustic signal database for a diurnal lineage of apose- 
matic and cryptic/non-defended taxa, the poison frogs. First, we showed 
that aposematic and non-aposematic species share similar extinction rates, 
and aposematic lineages diversify more and rarely revert to the non- 
aposematic phenotype. We then characterized mating calls based on 
morphological (spectral), behavioural/physiological (temporal) and 
environmental traits. Of these, only spectral and temporal features were 
associated with aposematism. We propose that with the evolution of anti- 
predator defences, reduced predation facilitated the diversification of vocal 
signals, which then became elaborated or showy via sexual selection. 
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1. Introduction 

Acoustic signals of courting males are often conspicuous to eavesdropping ene- 
mies, and these males face a trade-off between attracting partners and avoiding 
predators [1]. Diurnal signallers are especially vulnerable to predation owing to 
their increased visual and acoustic detectability [2]. These individuals may 
avoid predation with a Spectrum of strategies from crypsis to aposematism 
(the linking of a waming signal with a defensive strategy) [3]. Non-defended 
individuals usually signal from concealed places and rely on crypsis to avoid 
detection. If predation pressure is strong, individuals may reduce the number 
of signals, alter their Schedule or switch communication channels [4]. By con- 
trast, aposematic individuals (aposemes) broadcast warning signals that may 
be detected by predators [5]. Consequently, the evolution of aposematism, 
with the resulting reduction in predation risk, may enhance the diversification 
of mating signals and signallers. We tested this hypothesis by a phylogenetic 
analysis of the vocalizations of aposematic and cryptic species of the poison 
frog family Dendrobatidae. 

Aposematism is a complex phenotype that links conspicuous signals with 
defence (e.g. alkaloids). The association of these traits with other ecological 
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adaptations, including metabolic and diet specializations, 
characterizes the aposematic syndrome in dendrobatids [6]. 
Visual cues such as chromatic conspicuousness are often a 
warning signal to predators such as avian predators on 
poison frogs [7]. By contrast, acoustic cues wam predators 
such as bats, which avoid the ultrasound chirps from defended, 
noctumal moths [8]. Although both visual and acoustic signals 
may alert predators, their combined effects might increase the 
effectiveness of anti-predator defences [9]. 

Some predators locate their prey by exploiting components 
of the prey's signal. The efficacy of visual warning signals 
depends on light conditions and background [5]. Conspicuous 
prey might be effectively cryptic if the incident light is too low 

[10] . By contrast, sound-oriented predators and parasites 
recognize long-range signals such as mating calls and when 
in proximity switch to visual or Chemical identification [4]. 
For example, Corethrella midges locate male frogs by their 
calls and then switch to olfaction or other senses to locate the 
nostrils of the male, from which they obtain a biood meal 

[11] . Natural selection should favour individuals that avoid 
predator attacks by maximizing aposematic conspicuousness. 

Two general types of mechanisms may explain the origin of 
aposematism: first, predator-related mechanisms, including 
prey aggregation, dietary conservatism of predators and neo- 
phobia [5] and second, traits shaped primarily by natural 
selection, which can then be co-opted as sexual ornaments 
through sexual selection [12]. Aposemes might be better at 
attracting mates if predators associate their mating signals 
with unprofitability and thus avoid courting individuals. 
Thus, under reduced predator pressure, aposemes may 
evolve more easily detected mating signals via sexual selection. 

Poison frogs are a model clade for studies of the preda- 
tor -prey ecology of aposematism, which evolved at least 
four times in this group (figure 1). However, not all dendro- 
batids fit the stereotype of the brightly coloured, charismatic 
frog. Most are cryptically coloured and rely on camouflage. 
By contrast, the aposematic species are visually conspicuous 
and defended by skin alkaloids [13] that are distasteful and 
at times toxic to predators (e.g. birds, crabs and snakes). Den- 
drobatids are mostly diurnal and use visual and acoustic 
signals for intraspecific communication [14]. Their vocaliza- 
tions are innate and highly stereotyped, but have easily 
quantifiable variation among species. Within some species, 
females prefer males with greater calling performance [15,16]. 
Calling has a significant metabolic cost for male frogs [17] 
and aposematic dendrobatids have higher metabolic rates [6], 
but the relationship between acoustic signalling, aposematism 
and metabolic rates has not been explored in dendrobatids. 

Using phylogenetic methods, we tested whether acoustic 
mating signals diversified in association with the multiple 
origins of aposematism. Our results support this postulate 
and provide evidence that aposematism is associated with 
increased speciation rate in dendrobatids. We propose that 
the origin of aposematism, with the resulting reduction in 
predation risk, enabled mating calls of defended species to 
diversify via sexual selection. 


2. Material and methods 

(a) Acoustic data and perching behaviour 

We collected 16 657 advertisement calls from 172 species. All 
recordings were obtained from field collections and museum 


archives (electronic supplementary material. Dataset SI). The 
calls are characterized by single pulses with little frequency 
modulation (figure 1; electronic supplementary material, SI 
and Dataset SI). Field recordings were digitized with a sampling 
rate of 16 bits at a rate of 22 or 44 kHz and filtered for back- 
ground noise using a bandpass filter of 1-5 kHz. Spectrograms 
and power spectra were estimated using a Fast Fourier Trans- 
form (FFT) analysis using a Blackman window, 900 samples of 
overlap among subsequent FFTs, and 3 dB filter bandwidth of 
87.5 Hz. Homology of acoustic units was assessed following a 
physiological definition in which the call is the sound unit pro- 
duced by a cycle of trunk muscle contraction resulting in an 
expiratory event [18]. We used note-pulses, which have uniform 
temporal, spectral and taxon-specific features, as the homologous 
acoustic units. Temporal features and spectral properties were 
measured from oscillograms, spectrograms and power spectra 
(electronic supplementary material, figure S2 and Text). All 
acoustic variables were measured using RavenPro v. 1.4 [19]. 
We analysed 18 call variables measured from note-pulses (hom- 
ologous acoustic units) as well as temperature recorded at the 
calling site and body size (electronic supplementary material, 
figure S2 and table SI). For each call variable, the mean of the 
individuals was used for analysis. Finally, we also qualitatively 
described perch (calling) site as exposed or concealed based on 
published and direct observations for 83 species (electronic 
supplementary material. Dataset SI). 


(b) Alkaloid sequestration and conspicuousness 
variables 

We compiled skin alkaloid information of 97 taxa (electronic sup- 
plementary material, table SI). Species were characterized by 
their ability to sequester alkaloids as state 1 (able to sequester) 
or 0 (unable to sequester) [6]. A species was characterized as 
aposematic by the presence of defensive dietary alkaloids [6] 
and visual conspicuousness (electronic supplementary material, 
tables SI and S2). 

Chromatic contrast against a natural background is con- 
sidered a measurement of conspicuousness to predators [6]. 
Few dendrobatids have been assessed for conspicuousness 
using direct approaches such as total reflectance flux and 
models of predator perception [3,7]. Quantifying conspicuous- 
ness in life for 172 taxa using direct techniques was intractable. 
Therefore, we measured relative conspicuousness based on 
human perception of colour contrast against a leaf litter back- 
ground. Some authors [20] validly criticize the quantification of 
colouration based on human perception. These criticisms do 
not necessarily invalidate our analyses because most receivers 
include a mixture of trichromatic conspecifics and di-, tri- and 
tetrachromatic predators (electronic supplementary material, SI 
for discussion). All these receivers also have visual sensitivities 
that overlap with the human vision range (400-700 nm) and 
may not obtain information from the UV range [21]. 

As a proxy for direct approaches, we formulated a binary 
assessment of conspicuousness against a leaf litter background 
(electronic supplementary material, figure S3 and table S2). 
Colour descriptions of live male specimens were quantified by 
multiple independent human observers (X = 2.8 + 1.17 s.d. 
per species; electronic supplementary material, table S2). 
Using 11 frog skin segments, chromatic contrast (i.e. different 
from grey, brown and black) was scored as 1 (conspicuous) or 
0 (cryptic). The total contrast score (TCS or SS/), which 
ranged from 0 (no contrast) to 11 (maximum contrast), was 
determined by summing all 12 binary values. To account for 
inter-ob server variation, we used six cut-off values (2S; > 3, 
SS* > 4, . . . ,2S; > 8; electronic supplementary material, figure 
S3) of increasing colour contrast thresholds. These thresholds 
ranged from liberal (a species with XS* > 3 (TCS3) is categorized 
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Figure 1 . Poison frog chronogram showing the diversity of mating calls. Next to each tip (species) is the first 2 s of its spectrogram. Alkaloid sequestration is 
indicated by coloured dots: presence (orange), absence (brown) or unknown (grey). Red dots indicate species considered conspicuous based on the total contrast 
score 25/ > 5 (TCS5); black dots indicate cryptic species (25/ < 5). White dots indicate species calling from above-ground locations or in exposed habitats; purpie 
dots indicate calling from concealed locations. Black stars indicate origins of aposematism. Green dots indicate relative body size. Branch colours represent parsimony 
reconstructions of the temporal call component (PC2) using standardized units (bluish: slow rates; reddish: fast rates). 4//., Ållobates ; Ame., Ameerega ; Anom., 
Anomologlossus ; Col., Colostethus ; Rhe. f Rheobates ; Silverst., Silverstoneia. Support for each node is provided in electronic supplementary material, figure SI. 


as conspicuous) to strict (2 S/ > 8 (TCS8)). Thus, conspicuous- 
ness in each species was quantified under six binary variables 
TCS3-TCS8. For example, a given taxon with TCS7 is con- 
sidered cryptic (state 0) under the threshold 2 S/ > 8 (TCS8); 
but conspicuous (state 1) under the five thresholds 2 S/ > 3, 
2S/ > 4, . . . ,2 Si > 7 (i.e. TCS3-TCS7). This suite of variables 
assesses the robustness of our results over a range of receivers 
and light conditions. Finally, we also determined the ability of 
our conspicuousness variables to predict aposematism using a 
joint criterion for binary classifiers as described in the 
electronic supplementary material, SI. We emphasize that our 
conspicuousness assay is not an ideal substitute for modelling 
perception by predators. 


(c) Metabolic rate variables 

In electronic supplementary material, table SI, we compiled 
metabolic rate parameters of 54 species of poison frogs from a 
previous study [6] along with our corresponding acoustic data. 
The variables were (i) resting metabolic rate (RMR, oxygen con- 
sumption while resting or V0 2 rest ml h 1 ), (ii) active metabolic 
rate (AMR) after non-sustainable exercise (oxygen consumption 
after forced activity or V0 2 active mlh -1 ) and (iii) mean body 
mass to the nearest 0.01 g of all the individuals tested. From 
these raw data, mass-specific metabolic rates (AMR and RMR) 
were estimated by dividing the metabolic rates by the body 
mass of each individual. The average of all conspecific rates 
was used as the species mass-specific metabolic rate. Finally, 
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body mass and mass-corrected metabolic rates were transformed 
using natural logarithms to improve statistical distribution 
properties for the comparative procedures [22]. 


(d) Phylogenetic and comparative analyses 

The phylogeny was inferred from new and published molecular 
data: approximately 2400 bp 12S-16S rDNA mitochondrial 
genes (electronic supplementary material, table SI for GenBank 
numbers). The statistics of the molecular data matrix were as 
follows: (i) total sequence length (N = 172, X = 2330.13 + 
293.54 bp, missing cells are 13 910/414 692 or 3.35%); 
(ii) total sequence length per rRNA gene (12S: N = 172, 
X = 903.84 + 105.89 bp, missing cells are 3.44%; 16S: N= 172, 
X = 1358.88 + 175.42 bp, missing cells are 4.03%); and (iii) total 
sequence length per tRNA gene (f Val : N = 172, X = 67.41 + 
12.89 bp, missing cells are 7.65%). Tree estimation and nodal 
support were calculated under maximum-likelihood (ML) and 
Bayesian approaches using partitioned models. ML and Bayesian 
analyses gave similar tree topologies and the ML tree was used as a 
starting tree topology for the time-calibrated tree (electronic sup- 
plementary material, SI). The chronogram (figure 1; electronic 
supplementary material, SI) was determined using BE AST 
v. 1.5.3 [23] with five node-age constraints (electronic supplemen- 
tary material, figure SI). Based on these analyses, two taxonomic 
changes are made: Ameerega erythromos and Colostethus jacobuspe- 
tersi as part of Hyloxalus (i.e. H. erythromos and H. jacobuspetersi ; 
new combinations; see electronic supplementary material, SI for 
details). Chronogram tree file is deposited in the TreeBASE 
database under the accession number 16380. 

To quantify the relationship between components of apose- 
matism and call variables, we used diversification analyses, 
multivariate data exploration, tests of phylogenetic signal and 
models of trait evolution, bivariate phylogenetic correlations, 
exploratory factor analyses (phylogenetic principal component 
analysis, PPCA) and phylogenetic logistic regressions (PLRs). 
Using the binary variables alkaloid sequestration and conspicu- 
ousness, we estimated the rates of speciation (A 0 and Ai), 
extinction (/jl 0 and ^ i ) and transition between character States 
(qoi and q^o) using Binary State Speciation and Extinction 
(BiSSE) models [24]. Multivariate data explorations and variable 
reduction were used to narrow the dataset to 18 variables (169 
taxa) that loaded on three principal components. 

The PLRs were used to determine if call variables significantly 
predicted alkaloid sequestration and conspicuousness as depen- 
dent binary variables. Our predictors were two sets of continuous 
variables: the three PCs (principal components) derived from the 
PPCA of the call characters, and individual call variables, with 
body size and temperature as covariates. PLRs were performed 
with the PLogReg routine [25], which tests for phylogenetic signal 
while simultaneously performing the regressions. We applied the 
percentage increase in odds and the 'divide-by-4 rule' (i.e. y8/4 
where \ 3 is the regression coefficient) to determine significance of 
logistic regression coefficients [26]. See electronic supplementary 
material, SI for an example of y8/ 4 interpretation. Outliers were 
identified using standardized residuals with absolute values more 
than 3.0 and Cook's distance more than 1.0 as criteria. For pairwise 
correlations between the discrete dependent variables, we used 
Pagel's 1994 test for correlation of two binary characters [27]. Signifi- 
cance of all analyses was determined at a = 0.05, two-tailed 
distribution. 


3. Results 

We tested if acoustic courtship cues are associated with the 
origin of the aposematic phenotype, and if these pattems 
are related to species diversification. Specifically, (i) we 


determined whether conspicuousness and alkaloid sequestra- 
tion, the components of aposematism, are correlated and 
show phylogenetic signal; (ii) we measured species diversifi- 
cation by comparing extinction, speciation and character-state 
transition rates between aposematic and cryptic species; 
(iii) we used PPCA to describe the relationship between 
call variables and each aposematic component; and (iv) we 
used PLR to determine which call variables are associated 
with aposematism. 

In analysis (i), we used PageTs A [27] to test for phylogenetic 
signal (A > 0) in each of the two aposematic components scored 
(electronic supplementary material, figure S4 and table S3). Alka- 
loid sequestration showed phylogenetic signal (p < 0.001), 
as did conspicuousness variables TCS3-TCS6 (all p < 0.001). 
Alkaloid-bearing species were mostly conspicuous (all TCS vari- 
ables; all correlations significant at p < 0.032). Thus, aposematic 
species tend to be closely related and not randomly distributed 
across the phylogeny. 

In analysis (ii), we used BiSSE models [24] to estimate 
rates of speciation (A), extinction (pi) and transition between 
alternative character States (q). Given the concern about the 
interpretation of BiSSE models with less than 300 terminals 
[28], we assessed statistical power (electronic supplementary 
material, SI). As determined by simulations, our sample size 
(172 terminals) was adequate for analyses of speciation and 
transition rates, but less so for extinction rates. 

The speciation rate for conspicuous lineages showed a 
1.36- to 2.18-fold increase over that in cryptic lineages, 
except for TCS6 (p<0.05 for all; figure 2; electronic sup- 
plementary material, S4). However, for alkaloid-bearing 
lineages neither the speciation rate nor the extinction rate 
was different from that of non-defended lineages (p > 0.05 
for all). Given that the speciation rate for conspicuous 
lineages was higher, it is surprising that defended clades 
did not have a high speciation rate, because conspicuousness 
and sequestration are generally highly correlated. This result 
is perhaps explainable by the large amount of missing data 
for alkaloid sequestration (figure 1). However, it does not 
alter the general conclusion that aposematic clades have a 
higher speciation rate. The extinction rate for conspicuous 
lineages was not different from that for cryptic lineages 
(p > 0.05 for all) except for TCS3. 

We found that the rate of change from inability (state 0) to 
ability (state 1) to sequester alkaloids is 14 times higher than 
the reverse (q 01 = 0.014 versus q 10 < 0.001; p = 0.044; figure 2; 
electronic supplementary material, S4); essentially no rever- 
sals from the defended to the non-defended phenotype 
have taken place. Similar ly, the transition rate from cryptic 
to conspicuous States (q 01 ) and the reverse (q 10 ) shifted from 
a higher rate (q 01 > q 10 in TCS3 and TCS5; all p < 0.011) 
to a lower rate (q 01 < q 10 in variables TCS7 and TCS8; all 
p < 0.007). This is supported by two more pieces of evidence: 
TCS5 was significantly correlated with alkaloid sequestra- 
tion (p < 0.001), and TCS5 was also the best predictor of 
aposematism based only on conspicuousness (electronic sup- 
plementary material, figure S4). However, these results 
suggest that conclusions about rates of state change depend 
on how the binary state is defined, such as when large tip- 
ratios (more than 10 : 1) exist for change between States [28]. 
In summary, the extinction rates of aposematic and cryptic 
clades are not distinguishable, but aposematic lineages 
speciate more and are unlikely to revert to the cryptic/ 
non-defended phenotype. 
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Figure 2. Diversification analyses of aposematism in poison frogs. (a) Distribution of conspicuousness variables (TCS3-TCS8) derived from the total contrast score (TCS or 
25/) and alkaloid sequestration ability. PCI and PC2 represent call parameter space. Each species was coloured based on its conspicuousness and alkaloid sequestration 
ability. The Batesian mimic species (open cirdes) and the aposeme model (open squares) are also indicated. (b) Distributions of the rates of speciation- A, extinction — pc 
and transition between character States— c/. Alternative States are indicated by subscripts: O (cryptic/unable to sequester alkaloids) and 1 (conspicuous/able to sequester). 
The p-values refer to the probability of the nuli model of equal diversification rates (i.e. A 0 = A 1; /jlq = or p 10 = q o: ). 


In analysis (iii), PPCA recovered three components 
describing 77.88% of the variability (electronic supplemen- 
tary material, table S3): PCI (37.46% of the variance) 
described body size and spectral properties such as pitch, 
which are related to morphology [29]; PC2 (28.34%) 
described temporal properties (e.g. timing of pulsed notes), 
which are related to behaviour/physiology (e.g. contraction 
of trunk muscles for sound production) and PC3 (12.08%) 
reflected ambient temperature. All three PCs showed phylo- 
genetic signal (all A > O, all p < 0.001), indicating that call 
properties track phylogeny [18]. 

In analysis (iv), PLRs were used to determine if call 
components are associated with aposematism (figure 3; elec- 
tronic supplementary material, table S4). The predictor 
variables were the three PCs, and alkaloid sequestration 
and conspicuousness were discrete dependent variables. 
PCI (morphology) and PC2 (behaviour/physiology) pre- 
dicted alkaloid sequestration (all regression coefficients at 
p < 0.05). All else being equal (assumed hereafter), species 
with larger body sizes and low-pitch calls (PCI) have a 
20.2% ( + 6.18% s.e., p < 0.024) or a 1.2-fold increase in the 
probability of sequestering alkaloids per each standard unit 
increase of their PCI scores. Similarly, we found that species 
with short-duration pulses and/ or pulses at faster rates (PC2) 
have a 39.6% (±9.70% s.e., p < 0.018) or a 1.4-fold increase in 
their probability of sequestering alkaloids. The environ- 
mental component (PC3) was not a significant predictor 
(p>0.05). Therefore, species that call at low frequencies 
and /or emit vocalizations at faster rates are more likely to 
sequester alkaloids. 

Next, we investigated if PCI was associated with 
conspicuousness. PCI variables predicted conspicuousness 


only for variables TCS5 and TCS6 (p < 0.003 for both) 
but not the other variables. Species with larger body size 
or lower-pitched calls had approximately 1.2-fold increase 
(15.1 ± 5.03% s.e., p = 0.003 for TCS5 and 22.1 + 6.15% 
s.e., p < 0.001 for TCS6) in the probability of being conspicu- 
ous. Thus, these results support that body size and call pitch 
are associated with conspicuousness, and thus aposematism. 

These PLR results for PCI further support the associa- 
tion between body size, metabolic rates and aposematism 
as found previously [6]. However, is there a relationship 
between call characteristics and metabolic rates in poison 
frogs? We tested if mass-specific metabolic rates (AMR and 
RMR) are correlated with the acoustic trait PCs using a phy- 
logeny of 54 species having metabolic measurements [6] and 
acoustic data (electronic supplementary material, table SI). 
AMR and RMR were correlated with PCI (r AMR = 0.31 and 
Trmr = 0.43; both p < 0.001), reinforcing the connection of 
metabolic rate with spectral call properties (electronic sup- 
plementary material, tab le S3). Interestingly, metabolic rates 
were not correlated with PC2 or PC3 (both p > 0.05). Given 
that PC2 reflects performance-related variables such as 
pulse rate, our results reflect the greater stamina in body 
wall muscles (for calling) compared with limb muscles (for 
locomotion) as evidenced in hylid frogs [30,31]. 

The correlation of alkaloid sequestration with PC2 (tem- 
poral features) that was uncovered by our PLR analyses is 
novel. We further examined this by examining the individual 
variables that contribute to PC2. After controlling for body 
size and temperature, we found that pulse rate was the best 
predictor of alkaloid sequestration. An increase of one stan- 
dard unit resulted in a 36.2% ( + 9.55% s.e., p < 0.001) or a 
1.4-fold increase in the probability of the sequestration ability. 
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Figure 3. Phylogenetic logistic regression curves for alkaloid sequestration (top row) and conspicuousness (bottom row). The columns correspond to the principal 
components PCI: morphology (body size and spectral call parameters), PC2: behaviour and physiology (temporal call parameters) and PC3: environment. In the top 
row, solid lines are logistic regression curves for alkaloid sequestration, with 95% Cl (confidence interval; dashed lines) and (3 is the logistic regression coefficient. 
The bottom row shows logistic regression curves for each conspicuousness variable (TCS). Significant coefficients are marked with an asterisk (a = 0.05, two-tailed). 
Frequency histograms indicate the number of species along each principal component axis, in each category (able to sequester alkaloids or unable, conspicuous 
or cryptic). 


Shorter pulse duration, rise time and pulse intervals showed 
similar approximately 1.3-fold increases (p < 0.001 for all). 
Thus, individual temporal call variables (PC2) are sufficient 
to predict alkaloid sequestration. 

The individual variables of PC2 also predicted con- 
spicuousness (figure 3; electronic supplementary material, 
tab le S4). Species with short-duration pulses and /or faster 
pulse rates had a 12.0-27.1% (p < 0.05 for all conspicuousness 
variables) or approximately 1.2-fold increase in the probability 
of being conspicuously coloured. Individual temporal vari- 
ables also predicted conspicuousness. After controlling for 
body size and temperature, pulse rate had the highest predic- 
tive ability. An increase of one standard unit of pulse rate 
resulted in a 13.6-23.5% (all variables p < 0.05) or approxi- 
mately 1.2-fold increase in the probability of having the 
conspicuous phenotype. Similar results were found for other 
PC2-loading variables (electronic supplementary material, 
tab le S4). Thus, temporal features of mating calls are associated 
with conspicuous colouration. 

4. Discussion 

We propose two hypotheses for the connection between 
mating call characteristics and aposematism. Under the first, 
novel features of acoustic signals (shifts in spectral and tem- 
poral properties) evolved first as the primary warning 
component of aposematism and enhanced predator avoidance 
through experience with distasteful alkaloids. Visual conspicu- 
ousness (colouration) followed. However, we did not find 
phylogenetic or ecological support for this prediction. In a phy- 
logenetic context, the acquisition of novel acoustic signals 
should precede the appearance of the visual signal (i.e. conspic- 
uous colouration). However, we found the reverse patt em, in 
which changes in the vocalizations towards faster and lower- 


pitched calls occurred only after the appearance of conspicuous 
colouration (figure 1). 

Our first hypothesis also makes predictions about the 
ecology of the warning signal. Acoustic signals facilitate 
social interactions, but they also increase conspicuousness 
to predators. A predator might perceive the mating calls of 
defended individuals as a warning even in the absence of a 
visual warning cue. Many aposematic dendrobatids have 
conspicuous trills, chirps and buzz calls (electronic sup- 
plementary material. Dataset SI). Similar sounds such as 
the buzzing of bumblebees, the ratties of venomous snakes 
and the chirps of distasteful moths are warning signals [32]. 
Thus, we might expect aposematic dendrobatids to have 
calls with a broader frequency content Spectrum ('harsher') 
than those of non-aposematic species, which should have a 
narrower bandwidth or be pure tone ('sweeter'). Compari- 
sons between aggressive and attractive calls in birds and 
mammals [33] showed that these animals use harsh sounds 
when hostile and pure tone-like sounds when appeasing 
other individuals. Therefore, the calls of aposematic frogs 
might motivate avoidance by predators. However, we 
found no support for a correlation between broadband calls 
(measured by the bandwidth of the interquartile range) and 
either conspicuousness or alkaloid sequestration (both 
regression coefficients p > 0.05; electronic supplementary 
material, table S4). Thus, the first hypothesis, that the 
mating calls of aposematic dendrobatids are a primary 
warning signal, is not supported. 

Our second hypothesis proposes that the evolution of 
aposematism, with the resulting reduction in predation risk, 
may have freed mating calls to evolve showy characteristics 
via sexual selection. This hypothesis is favoured by both phy- 
logenetic and ecological evidence. In the evolutionary context, 
aposematism preceded the evolution of the diversification 
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mating calls (figure 1). Consequently, any change in call traits 
occurred after aposematism was already a functional anti- 
predator strategy. In the ecological context, several aspects of 
the ecology of dendrobatid aposemes support our postulate. 
The diurnal activity of almost all dendrobatids includes multi- 
modal courtships that integrate acoustic and visual displays 
[14]. Visually oriented predators learn to avoid aposemes, 
which then undergo ecological release and evolve existing 
calls into more conspicuous displays that attract females. 
Natural history information and our direct observations 
indicate that most non-aposemes vocalize from concealed, 
ground-level sites, whereas aposemes vocalize from exposed, 
above-ground sites. We found that exposed callers are more 
likely to sequester alkaloids than concealed callers (Pagel's 
1994 test; electronic supplementary material, table S3, N = 
83, p = 0.015). Similarly, we found that exposed callers are 
more likely to be conspicuous than cryptic, concealed callers 
(electronic supplementary material, table S3, TCS3-TCS7, 
N = 146, all p < 0.020). 

These results suggest that calling behaviour of aposemes 
might result in higher efficiency of acoustic transmission. In 
tropical frogs, males that broadcast above the ground in 
open environments and call with lower pitch suffer less 
degradation over distance [34]. These observations agree 
with our findings that aposemes call at a lower pitch and 
also from higher perching sites than do cryptic, undefended 
frogs. Consequently, the calls of aposematic males should 
be more salient to potential mates, and thus females would 
prefer males that call from exposed perches. Therefore, this 
second hypothesis is better supported and suggests that 
aposematism releases predator-induced constraint on the 
evolution of the signal and the behaviour of the signaller. 

Three possible mechanisms explain how females might 
evaluate aposematism in mate choice. First, females select 
males that provide them direct benefits during courtship, 
such as effective predator deterrence. Second, males are 
selected if they provide indirect benefits to offspring such 
as increased viability and attractiveness derived from traits 
associated with aposematism. Third, males are selected due 
to sensory exploitation of aposematic signals by females. 
No experimental evidence exists for these mechanisms, but 
some inferences might be derived from comparison of field 
and laboratory studies. For example, males might gain 
mating advantages from female preference for conspicuous 
colouration, which originally evolved in a non-sexual context 
to warn predators. However, under certain laboratory con- 
ditions, females of Dendrobates pumilio prefer brighter and 
possibly more toxic males [35], while in natural conditions 
females prefer the closest calling male [36]. 

The exploitation of aposematic signals by mimics might 
also confer reproductive advantages to them and drive the 
evolution of their mating signals. Allob ates zaparo , a Batesian 
(e.g. non-defended) mimic, derives protection from matching 
the reddish dorsal colour of its aposematic models Ameerega 
bilinguis and Ameerega parvula [37]. Interestingly, we found 
that A. zaparo does not acoustically resemble either model 
closely (figure 2a; electronic supplementary material, 
figure S4A), but rather is acoustically doser in call parameter 
space to its non-aposematic relatives within the Allobates 
femoralis complex (Euclidean distance: 0.809 + 0.253 s.d.). 
However, calls of A. zaparo are within the range of variation 
of most aposemes; the PC scores of A. zaparo were not outliers 
in the PLR analysis (all Cook's distances < 0.042). This 


suggests that A. zaparo might benefit from the evolution 
of showy mating calls as an advantage of mimicking 
aposematic colouration and the resulting reduced predation. 

Acoustic signalling increases detectability, suggesting that 
predation pressure is greater on males [12] than on females, 
which have limited acoustic signalling (e.g. release calls). 
However, given that sexual dichromatism is not common in 
aposematic dendrobatids and the shared genetic architecture 
is present in males and females, predation pressure rather 
than mate choice is probably the main driver of visual con- 
spicuousness. In one species (D. pumilio ), females have 
larger quantities of alkaloids than males [38]. Hence, why 
may we expect aposematic females to have better defences 
than males? Some aspects of female ecology suggest a 
broader pattern of social selection. First, females are exposed 
to predators as they search and compete for the best partner 
and share with males the cost of parental care [39]. Second, 
the cost of aposematism might be higher if females provide 
defensive alkaloids directly to their offspring [40]. Third, 
dendrobatids engage in prolonged and interactive courtships 
where females search for and move towards vocalizing males 
[14]. This behaviour might increase females' exposure to pre- 
dators compared to males', but this hypothesis needs further 
experimental evidence. Finally, females vigorously defend 
general-purpose territories, attack intruders, and may use 
aposematism to protect ecological resources [41]. Further 
research on the behavioural ecology of females rather than 
only phylogenetic comparative analyses (i.e. correlation 
does not necessarily means causation) will illuminate the 
sex-specific benefits and costs of aposematism. 

In sum, our analyses revealed a strong association 
between acoustic signalling and aposematism in poison 
frogs. We provide two general conclusions. First, aposematic 
lineages speciate at a higher rate than non-aposematic 
lineages and rarely revert to a cryptic phenotype. A high spe- 
ciation rate was found in the aposematic clades Dendrobates 
and Ameerega [42], which evolved aposematism indepen- 
dently. This is consistent with the suggestion that female 
preference for males with brighter visual signals in D. pumilio 
[35] might drive species diversification. However, this 
assumes intra-populational variation in brightness, which is 
not demonstrated in most species. Nonetheless, our results 
might suggest that the higher speciation rates in aposematic 
lineages are due to increased sexual selection once predator 
pressure is diminished. 

Second, lower pitch (spectral properties) and faster rates in 
pulse-related call variables (temporal properties), as well as 
calling from exposed perches, are associated with multiple ori- 
gins of aposematism. Consequently, our results show that 
dendrobatid aposemes shifted to a different type of call com- 
pared to non-aposemes. This pattern is convergent among 
aposematic clades; thus dendrobatid calls have become more 
diverse within Dendrobatidae. We proposed two hypotheses 
for the evolutionary shift of acoustic signals and behaviour of 
aposemes. First, this shift was due to the superposition of a 
warning function onto the pre-existing mate attraction func- 
tion, which when combined with existing alkaloid defence, 
enhanced protection against predators. This suggests that 
receivers need only hear the frog to know if it is aposematic 
and hence unpalatable. However, our analyses do not support 
this postulate and it remains to be tested whether calls alone of 
aposemes deter predators. Second, this evolutionary shift 
resulted from pre-existing aposematism (alkaloid defence 


rspb.royalsocietypublishing.org Proc. R. Soc. B 281 : 20141761 


Downloaded from rspb.royalsocietypublishing.org on October 15, 2014 


and visual conspicuousness), which relaxed predation 
pressure. The combination of low-pitch and exposed calling 
sites enhanced call transmission and mate attraction. This 
last hypothesis is supported by our data and suggests that, in 
general, predation exerts a strong constraint on sexual 
selection; aposematism removes such constraint. 

Although further experiments will determine to what 
extent specific call traits associated with aposematism affect 
mate choice, our data offer evidence of the effect of preda- 
tion on the evolution of signal complexity. While a similar 
pattem of predators driving phenotypic diversification 
has been demonstrated in other taxa such as Trinidadian 
guppies (. Poecilia reticulata ) and Gambusia mosquitofish 
[43,44]. Our results pro vide evidence for the interaction of 
signals across modalities (visual and acoustic) and across con- 
texts (signalling to predators and to mates) to affect the 
evolution of communication, as well as the formation of 
new species in aposematic clades. Only studies of phenotypic 
predictors of reproductive success would be required to fully 
demonstrate that sexual selection via female choice is respon- 
sible for the further diversification of elaborate vocal signals 
in dendrobatid frogs. 
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